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The  infrared  spectroscopy  and  rovibrational  dynamics  of  HC1  and  DC1  dopants  in  solid 
parahydrogen  (pH2)  were  investigated  using  high-resolution  spectroscopic  methods.  The 
absorption  spectra  of  HC1  and  DCI .monomers  in  solid  p>H2  closely  resemble  the  corresponding 
low  temperature  gas  phase- spectra,  indicating  that  the  gas  phase  vibrational  and  rotational 
quantum  numbers  of  the  impurity  are  conserved  within  the  pH2  solid.  Small  deviations  from  gas 
phase  behavior  are  observed,  however,  such  as  a  reduced  HC1  rotational  energy  level  spacing  and 
splitting  of  the  fivefold  orientational  degeneracy  of  HC1  rotational  states  with  J- 2.  In  addition, 
the  pure  vibrational  Qi(0)  (v=l<-0,  J=0*-0)  H2  transition,  which  is  infrared  inactive  in  pure  solid 
pH2,  is  detected  in  the  HC1  doped  sample.'  This  transition  is  induced  in  pH2  molecules  by 
neighboring  HC1  molecules  through  a  weak  .‘‘overlap  induction”,  mechanism  that  is  the  only 
induction  mechanism  operative  for  J=0  impurities  in  pH2.  Rovibrational  transitions  are  also 
detected  near  the  induced  QfO)  H2  absorption;  these  are  attributed  to  cooperative  transitions 
involving  single  photon  excitation  of  pH2-HCl  pairs.  Detailed  isotopic  analysis  reveals  that 
these  cooperative  transitions  involve  pure  vibrational  excitation  of  the  pH2  and  pure  rotational- 
excitation  of  the  HC1.  Two-molecule  transitions  have  long  been  studied  for  isotopic  and 
rotational  dopants  (e.g:,  D2,  HD,  orthohydrogen)  in  solid  pH2,  but  this  is  the.  first  time  such 
cooperative  transitions  have  been  attributed  to  a  chemical  impurity  in  pH2  matrices.^ 
Corresponding  author:  danderso@uwyo.edu 

bPresent  address:  KLA-Tencor  Corp.,  1  Technology  Drive,  Milpitas,  CA  95035 


I.  INTRODUCTION 

The  infrared  (IR)  absorption  spectroscopy  of  hydrogen  halide  dopants  HX  (X  =  F,  Cl,  Br) 

-  in- rare  gas  (Rg)  solids  has  been  studied  extensively,  beginning  in  the  1960s.1  These  studies, 
revealed  that  HX  dopants  in  rare  gas  solids  exhibit  R( 0),  i?(l),  and  PCI)  rovibrational  transitions, 
as  would  be  expected  for  HX  molecules  at  temperatures  of  10  to  20  K;  this  indicates  that  the 
rotational  motion  of  HX  molecules  embedded  in  rare  gas  matrices  is  only  weakly  hindered.  For 
HC1  dopants,  this  interpretation  was  later  confirmed  by  observations  of  both  the  HC1  .7=1  «-0 
pure  rotational  absorption  transition2  and  the  rovibrational  Q  branch  in  the  HC1  Raman. 
spectrum.3  There  are  two  main  differences  between  the  rovibrational  spectrum  of  gas  phase  HC1 
and  the  spectra  exhibited  by  HC1  dopants  in  rare  gas  solids,  differences  which  arise  from  the  Rg- 
HC1  intermolecular  potential:  a  slight  redshift  in  the  vibrational  band  origin  and  a  small 
reduction  in  the  rotational  spacing.  Analysis  of  the  spectral  perturbations  induced  by  the 
-cryogenic  matrix  can  therefore  provide  information  on  the  rovibrational  dynamics  of  HC1 
dopants  in  the  matrix  environment  and  the  potential  surface  felt  by  these  dopants  in  this 
environment.  -v  ...  ...  ...  ■  .... 

For  HC1  dopants  located  in  single  substitution  sites  of  a  face  centered  cubic  (fee)  rare  gas 
crystal,  the  octahedral  site  symmetry  effectively  eliminates  many  of  the  anisotropic  terms  in  the 
Rg-HCl  potential;  hence  a  dopant  HC1  molecule  with  its  center  of  mass  held  fixed  at  the  center 
of  the  substitution  site  feels  a  nearly  isotropic  orientational  potential  surface.  However,  if  the 
rare  gas  matrix  relaxes  around  the  HC1  dopant,  lowering  the  octahedral  symmetry  of  the 
substitution  site,  or  if  the  HC1  moves  away  from  the  center  of  die  cavity,  the  potential  anisotropy 
increases,  leading  to  energetic  barriers- to  the  rotation  of  the  HC1  molecule  in  the  matrix.  This 
perturbs  the  rotational  motion  of  the  HC1,  and  can  both  shift  the  energy  of  excited  HC1  rotational 
states  with  J>0  and  split  the  (2J+l)-fold  orientational  degeneracy  of  these  states.  This 
perturbation  of  the  HC1  rotational  quantum  states  by  angular  anisotropy  is  historically  called  a 
crystal  field  effect.4'5 

Another  perturbation  of  the  dopant  rotational  states  results  from  the  coupling  between  the  > 

dopanf^translational  and  rotational. piotion.6'8  As  the  HC1  rotates,  its  minimum  energy  center- 

\yijp  .  »  v  '  "  ;  ■  :  '  ,,a. 

of-mass  position  in  the  single  substitution  site  varies,  leading  to  a  coupling  between  the 

molecule’s  translational  and  rotational  degrees  of  freedom  as  it  both  rotates  and  rattles  around  in 
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the  rare  gas  cavity.  This  rotation-translation  coupling  can  also  lead  to  deviation's  of  the 
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condensed  phase  spectrum  from  the  gas  phase  spectrum.  ’ 

.With  these  studies  in  mind,  we  decided  to  investigate  the  infrared  spectra  of  HC1  and  DC1 
in  cryogenic  solid  hydrogen  matrices  to  study  the  rovibrational  dynamics  of  these  species  in  this 
environment.  While  even  neon  rare  gas  solids  are  fairly  well  described  as  classical  crystals,  the 
light  mass  of  H2  and  weakly  attractive  H2-H2  intermodular  forces  conspire  to  make  solid 
hydrogen  a  quantum  solid,9'12  with  physical  properties  that  are  dominated  by  large  amplitude 
zero-point  motions  of  "the  individual  H2  molecules  about  their  nominal  lattice  positions.  In 
addition,  because  the  H2  rotational  constant  is  large  and  the  anisotropic  components  of  the  H2-H2 
interaction  potential  are  small,  individual  H2  molecules  retain  their  gas  phase  vibrational  and 
rotational  quantum  numbers  in  the  solid.  The  physical  properties  of  these  quantum  solids  make 
them' ideal -“solvents”  in  which  to  study  the  high-resolution  spectroscopy  and  rovibrational 
dynamics  of  chemical  impurities  in  a  condensed  phase.  Accordingly,  solid  molecular  hydrogen 
has  been  the  subject  of  numerous  investigations,13'28  with  both  pure  and  doped  systems  receiving 
considerable  attention. 


Molecular  hydrogen  exists  in  two  distinct  nuclear  spin  modifications:  orthohydrogen,  in 
which  the  nuclear  spins  are  coupled  together  to  form  a  triplet  spin  state,  and  parahydrogen,  in 
which  the  nuclear  spins  are  coupled  together  to  form  a  singlet  spin  state.  Orthohydrogen  (0H2) 
^molecules  have  odd  rotational  quantum  number^J^vhile  parahydrogen  (phk)  molecules  have 
even  values.  A  pure  orthohydrogen  solid  is  made  up  of  a  collection  of  J=1  molecule^whidy^ 
due  to  the  threefold  Mj  orientational  degeneracy  and  the  permanent  quadrupole  moment  of  the 
J=  1  stat^an  form  orientationally  ordered.phases  9,10  By  comparison,  solid  parahydrogen  is  made 
up  of  J=0  molecules,  which  are  individually  spherically  symmetric  objects  lacking  any 
permanent  electrostatic  moments.  Thus  solid  parahydrogen  is  a  collection  of  spherical  objects 
held  together  only  by  isotropic  intermolecular  forces,  much  like  a  rare  gas  solid.  However,  an 
important  distinction  between  solid  pH2  and  rare  gas  matrices  is  that  the  H2  molecules  that  make 
up  Solid  pH2  have  both  rotational  and  vibrational,  degrees  of  freedom.  In  this  study,  we  exploit 
this  distinction  to  probe  the  microscopic  dynamics  of  both  the  pH2  host  and  the  HC1  or  DC1 
dopant,  dynamics  which  are  mediated  by  the  H2— HC1  intermolecular  potential. 

Solid  pH2  is  a  hexagonal  close  packed  (hep)  crystal10  in  which  nearest'  neighbor  H2 
molecules  are  separated  by  3.79  A.  The  equilibrium  distance  of  the  pH2-HCl  intermolecular 
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potential  surface  is  3.88  A  (see  below).  Consequently  we  expect  that  the  HC1  dopant  should 
replace  a  single  pH2  molecule  in  the  hep  crystal.  Due  to  the  high  symmetry  of  the  single 
substitution,  site,  in  which  the  spatial  distribution  of  the  pHi  molecules  around  the  HC1  is  very 
nearly  spherical,  the  anisotropic  orientational  potential  felt  by  the  HC1  due  to  its  twelve  nearest 
neighbor  pH2  molecules  should  be  exceedingly  small.  Hence  the  HC1  should  freely  rotate  with 
its  center-of-mass  located  (on  average)  at  the  center  of  the  site,  as  is  the' case  when  HC1  is  doped 

-.19' 

into  rare  gas  solids.  ’ 

In  this  paper  we  'demonstrate  that  HC1  and  DC1  impurities  do  in  fact  rotate  nearly  freely 
in  solid  pH2  matrices,  so  that  the  rotational  quantum  number  of  the  impurities  remains  good  in 
the  matrix  environment.  Thus  the  infrared  absorption  transitions  of  HC1  and  DC1  dopants  in  pH2 
.  can  be.labeled  using. gas  phase  notation.  We  adopt  the  standard  Lp (V ,’)  notation  for  labeling  (v’, 
7)  <—  (v”=0,  7”)  rovibrational  absorptions,  in  which  L  is  given  by  O,  P,  Q,  R,  or  S  for  AJ  =  J  - 
J”  —  -2,  -1,  0,  +1,  or  +2,  respectively.  We  find  that  7= 0  HC1  and  DC1  dopants  also  induce 
infrared  activity  in  neighboring  pH2  molecules  through  a  weak  overlap  induction  mechanism. 
Finally,  we  report  the  first  observations  of  cooperative  transitions  in  chemically  doped  pH2 
matrices,  in  which  a  single  photon  excites  the  rotational  motion  of  a  dopant  molecule  and  the 
vibrational  coordinate  of  an  adjacent  pH2  molecule.  These  two-molecule  cooperative  transitions 
are  the  condensed  phase  analogues  of  combination  bands  in  van  der  Waals  clusters. 

The  rest  of  the  paper  is  organized  as  follows.  In  Section  II  we  present  the  experimental 
details  germane  to  this  paper.  In  Section  III  we  present,  assign  and  partially  analyze  the 
spectroscopic  data  related  to  the  direct  infrared  excitation  of  the  HC1  and  DC1  impurities.  In 
Section  IV  we  present  and  analyze  our  observations  related  to  rovibrational  transitions  Of  pH2 
molecules  in  the  matrix,  including  the  pH2-HCl  cooperative  transitions.  In  Section  V  we  discuss 
the  rovibrational  dynamics  of  HC1  dopants  in  a  pure  v=0  pH2  solid  and  the  rotational  dynamics  of 
HC1  dopants  with  either  a  nearest  neighbor  v=l  pH2  molecule  or  a  nearest  neighbor  v=l ,  7=2  H2 
molecule.  Finally,  in  Section  VI  we  formulate  our  conclusions  and  discuss  future  research 
directions. 
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II.  EXPERIMENT 

The  experimental  apparatus  and  sample  preparation  techniques  have  been  described  in 
....  detail  elsewhere;26,27  here,  we  .will  give  only  a  brief  swmmary^  highlighting,  the. .-specific  aspects 
relevant  to  the  present  experiments.  Doped  pHb  solids  are  prepared  by  rapid  vapor  codeposition 
of  precooled  pEb  gas  and  room  temperature  HC1  (or  DC1)  gas  onto  a  BaF2  substrate,  which  is 
cooled  to  T  «  2  K  in  a  liquid  helium  (lHe)  bath  cryostat.  The  precooled  pH2  is  prepared  in  an 
ortho/para  H2  converter  operated  at  15  K  with  11H2  (normal  hydrogen)  inlet  rates  «  200  mmol/h, 
yielding  a  pH2  flow  that  impinges  upon  the  substrate  at  a  45°  angle.  The  flow  rate  of  the  HC1  gas 
is  adjusted  to  produce  samples  with  10-300  parts  per  million  (ppm)  of  HO  in  the  pH2  host 
matrix.  We  prepared  both  samples  doped  with  natural  abundance  HC1  or  DO  and  samples 
.  doped' with  isotopically  enriched  H35C1.  H37C1,  D35C1,  or  D37C1,  prepared  by  reaction  of 
—  isotopically  enriched  NaC^Pwith  concentrated  H2SO4  and  D2SO4.  The  HC1  and  DC1  doped  pH2 
samples  typically  contain  residual  0H2  (J=l)  impurities  at  concentrations  of  roughly  100  ppm. 
,'toa  S.K2  ^  However,  as  we  describe  below,  some  pH2  samples  were  intentionally  doped  with  higher 
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concentrations  of  0H2  by  operating  the  ortho/para  converter  at  elevated  temperatures. 

; .  The  infrared  spectra  of  the  doped  pH2  solids  were  recorded  at  resolutions  ranging  from 
vf  0.005  to  0.01  cm’1"  the  main  optical  axis  is- parallel  to  the  substrate  normal.  The  Fourier 

'U '  transform  (FT)  IR  spectrometer  (Broker  IFS  120HR)  is  equipped  with  a  glowbar  source,  a  KBr 

beam  splitter,  and  an  InSb  detector  (1800—9000  cm  ^ ).  To  accommodate  the  IR  spectroscopy,  the 
lHe  cryostat  resides  inside  a  0.5  m3  polycarbonate  box  purged  with  a  constant  flow  of  dry  N2  gas. 
As-deposited  spectra  are  -recorded  immediately  after  deposition.  Annealed  spectra  are  recorded 
after  raising  the  temperature  of  the  solid  to4  .-8  K  for  a  specified  period  of time. 

In  this  study,  we  depart  from  our  normal  convention  of  reporting  concentrations 
determined  via  the  application  of  Beer’s  Law  to  the  observed  IR  absorption  spectra.  We  lack  a 
single  spectroscopic  feature  [e  g.,  the  HC1  or  DC1  monomer  /?  1  (0)  line]  which  remains  observable 
over  the  wide  range  of  dopant  concentrations  employed.  Instead,  we  report  the  concentration 
■  calculated  as  the  quotient  of  the  inlet  quantities  of  dopant  ,  and  pH2.  We  estimate  that  these 
reported  concentrations  are  within  ±  40  %  of  the  true  dopant  concentrations. 


HI.  RESULTS  AND  ANALYSIS:  TRANSITIONS  OF  HCI  AND  DCI  DOPANTS 

We  begin  by  extracting  the  J=  0  H2  component  from  the  full  4D  H2-HClintermolecular 
potential  energy  surface  (PES)29  to  'understand  the  nature  of  the  HCI  substitution  site  in  pH2. 
matrices.  Figure  1  shows  a  contour  plot  of  the  resulting  pH2-HCl  intermolecular  potential  in  a 
molecule  fixed  Cartesian  coordinate  system.  The  origin  of  the  coordinate  system  is  the  HCI 
center  of  mass  and  the W-axis  coincides  with  the  HCI  bond  axis.  The  pair  potential  has  a  global 
minimum  Vmm=  -88.97  cm*1  at  5=3.88  A  for  the  linear  pH2-HCl  configuration,  and  a  local 
minimum  V=  -65.46  cm'1  at  R= 3.52  A  for  the  linear  pH2-ClH  configuration.,  Because  the 
nearest-neighbor  spacing10  in  hep  crystalline  pH2  is  3.79  A,  we  conclude  that  HCI  should  easily 
be  accommodated  in  a  single  substitution  site.  Given  the  isotropic  nature  of  a  single  substitution 
site,  the  modest  anisotropy  in  the  pH2-HCl  PES  (see  Fig.  1),  and  the  large  rotational 
constants30’31  of  HCI  (5=10.44  cm'1)  and  DCI  (5=5.39  cm'1),  the  total  angular  momentum 
quantum  number  J  associated  with  HCI  and  DCI  dopants  in  pH2  matrices  should  remain  a  good 
quantum  number. 

A.  Survey  scans  of  HCI  and  DCI  doped  pH2  samples 

-  Survey  scans  of  the  HCI  and  DCI  doped  samples  -in  the  vicinity-of  dopant  vibrational 
fundamental  are  shown  in  Fig. '2.  As  expected,  the  as-deposited  spectrum  of  HCl  in  solid  pH2  is 
dominated  by  the  5i(0)  absorption  feature  of  the  H35C1  and  H37C1  dopants  present  in  natural 
abundance  in  the  HCI  sample.  The  isotopic  label  associated  with  each  peak  is  easily  identified 
from  the  characteristic  3:1  isotopic  abundance- ratio  and  the  2.1  cm'1  frequency  shift30  between 
the  H35C1  and  H37C1  fundamental  features;  these  assignments  are  unambiguously  confirmed 
using  isotopically  enriched  samples  (see  below).  -  . 

Upon  annealing  the  sample  at  5=4.8  K  for  short  periods  of  time,  additional  sharp  features 
grow  in  irreversibly  at  energies  below  that  of  the  HCI  monomer  transition.  -These  sharp  features 
are  due  to  hydrogen  bonded  (HC1)„  clusters  that  agglomerate  in  the  pH2  solid  at  higher 
temperatures  following  thermally-activated  HCI  diffusion.  Furthermore,  due  to  the  characteristic 
shift  of  the  HCI  stretching  frequency  upon  hydrogen  bond  formation,  the  different  cluster  species 
(n=2,  3,  etc.)  absorb  at  different  frequencies  and  are  well  resolved  in  the  infrared  spectrum.  Note 
that  in  the  as-deposited  sample  the  dominant  species  is  the  HCI  monomer,  attesting  to  the 
isolation  efficiency  of  the  rapid  vapor  deposition  method.26,27  The  cluster  features  are  of  interest 
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due  to  .  the  pervasive  role  that  hydrogen  bonding  interactions  play  in  a  variety  of 
environments,32,33  but  further  discussion  of  these  transitions  is  beyond  the  scope  of  the  present 

paper  and  will  be  presented  separately  in  a  future  publication/4  ... 

'  A  similar  spectrum  is  recorded  for  DCI  doped  samples  in  the  vicinity  of  the  DCI 
fundamental.  Again,  the  spectrum  of  the  as-deposited  samples  is  dominated  by  the  DC1  R\(0) 
monomer  transition,  which  is  actually  off  scale  in  Fig.  2.  The  D3dC1  and  D37C1  isotopomers  are 
present  in  the  3:1  natural  isotopic  abundance  ratio  and  give  rise  to  two  features  separated  by  the 
3.0  cm'1  frequency  difference31  characteristic  of  the  fundamental  DC1  stretch.  Upon  annealing 
the  sample,  sharp  transitions  grow  in  irreversibly  with  comparable  intensity  patterns  and  spacing 
to  those  measured  in  the  HC1  doped  samples.  By  analogy,  these  additional  features  are  ascribed 
to  (DC1)„  clusters  that  form  in  the  pH2  solid  at  elevated  temperatures.  . 

For  both  HC1  and  DC1  dopants,  some  of  the  sharp  features  attributed  to  the  clusters 
appear  very  close  to  the  vibrational  band  origin  of  the  corresponding  monomer.  These  cluster 
.features  are  presumably  due  to  excitation  of  the  “dangling”  HC1  or  DC1  bonds  of  the  various 
cluster  species.  These  so-called  “free”  or  “acceptor”  HC1  and  DC1  molecules  are  not  shifted 
strongly  by  the  hydrogen  bond  interaction  and  absorb  near  the  monomer  origin.32,33  Regardless 
of  the  specifics  of  the*  cluster  assignment  for  these  features,  it  is  clear  that  they  can  (m  principl^ 
overlap  with  or  obscure  the  monomer  ro vibrational  transitions.  However,  comparison  of  the  as- 
deposited  and  annealed  samples  can  facilitate  tentative  assignments  of  features  to  isolated 
dopants  or  to  clusters;  these  assignments  can  be  tested,  via  further  analysis.  For  example,  in  the 
next  section  we  show  that  a  small  peak  at  roughly  2072.2  cm'1  that  is  just  discemable  in  the  as- 
deposited  DCI  sample  arises  from  the  £i(l)  transition  of  isolated  DCI  dopants. 


B.  DCl/pH2  rovibrational  assignment 

Assignment  of  the  DCI  -monomer  features  is  done  by  inspection  based  on  a  model  of 
nearly  free  rotation  of  DCI  dopants  in  the  pH2  matrix.  Based  on  the  gas  phase  rotational 
constant31  of  DCI  (£=5.39  cm'1  for  v=0),  at  the  highest  temperature  (£=4.8  K)  accessed  in  these 


experiment^approximately  89%  of  the  DCI  molecules  are  in  the  ground  J- 0  rotational  state, 
with  11%  in  J=\  and  less  than  0.03%  in  J=2.  Hence  the  rovibrational  spectrum  of  DCI  in  pH2 


should  be  dominated  by  a  single  i?i(0)  transition  at  low  temperature  with  i?i(  1)  and  P i(l) 


transitions'  visible  in  the  higher  temperature  scans. 


'  Otf. 


Infrared  spectra  of  isotopically  enriched  D35Cl/pH2  and  D37Cl/pH2  samples  are  shown  in 
Figs.  3  and  Respectively.  The  observed  transition  frequencies  and  linewidths  (full  width  at  half 
maximum  or  FWHM)  are  summarized  in  Table  I  for  both  isotopomers.  Trace  (a)  is  the  as- 
deposited  D35Cl/pH2  sample  in  which  contributions  from  (D35C1)„  cluster  absorptions  are 
minimal.  The  Pi(l)  transition  is  evident  in  this  spectrum  approximately  48  or  20  cm'1  lower  in 
energy  than  the  strong  i?i(0)  peak,  which  is  off  scale  in  the  .spectrum.  Traces  (b)  through  (d) 
show  spectra  recorded  at  T=4.8,  2.4,  and  3.6  K,  respectively,  demonstrating  the  reversible 
temperature  dependence  of -this  monomer  feature.  The  features  near  2078  cm'1  assigned  to 
(DC1)„  clusters  irreversibly  increase  in  intensity  during  the  initial  annealing.  The  “hot”  ?i(l)  and 
i?i(l)  rovibrational  transitions  at  roughly  2072.2  cm'1  and.  2102.6  cm'1,  respectively,  can  be 
assigned  based  on.  their  reversible  temperature  dependence  and  -4 B  and  +2 B  spacing  of  these 
features  from  R\(0). 

The  8i(0)  and  Pi(l)  transitions  display  no  visible  (<0.1  cm'1)  fine  structure  an^hu^hese 
two  lines  are  used  to  calculate  an  effective  rotational  constant  5Avg,  where  5Avg  =  (8y=o  + 
5v=i)/2.  This  leads  to  a  value  of  8AVg  =  5.1  cm'1  (experimentally  indistinguishable,  for  both 
D35C1  and  D37C1)  which  is  close  to  but  slightly  smaller  than  the  gas  phase  values31  of  5.34  cm'1 
for  D35C1  and  5.32  cm'1  for  D37C1.  The  slight  reduction  in  B  in  the  matrix  environment  arises 
from  crystal  field  effects,  as  explained  earlier.  Based  on  this  value  of  5AVg,  the  DC1  i?i(l) 
transition  should  be  observed  at  2102.8  cm'1.  Consistent  with  this  prediction,  a  “hot” 
rovibrational  transition  is  -detected  in  tins  region  that  is  split  into  (at  least)  three  components. 
The  average  frequency- of  the  three  resolvable  features  is  2102.6  cm'1  for  D35C1,  which  is  very 
close  to  the  predicted  £i(l)  "transition  frequency.  The  origin  of  the  splitting  of  the-  8;(1) 
transition  will  be  discussed  below.  With  this  rovibrational  assignment,  the  vibrational  band 
origin  for  D35C1  and  D37C1  are  calculated  to  be  v0  =  lA[Ri(0)  +  Pi(l)j  =  2082.4  cm'1  and  2079.4 
cm'1,  respectively,  each  redshifted  by  8.7  cm'1  from  the  corresponding  gas  phase  value.31 

To  further  verify  that  the  Pi(l)  and  i?i(l)  transitions  originate  from  a  thermally  populated 
.  J=1  DC1  state  at  approximately  Emt  =  10.2  cm'1,  the  temperature  dependence  of  the  integrated 
absorbance  of  these  features  was  fit  to  a  standard  Boltzmann  equation, 


J")  =  Nj~{J"+  J'+ 1)  exp(-2 B/kT)/Qmt(B,  T) 


where  Nj*  and  B  are  determined  via  a  nonlinear  least  squares  fit  to  the  experimental  data  shown 
in  Figs.. 3  and  4.  0:ot(B,  T)  is  the  rotational  partition  function  (calculated  for  a  rigid  rotor)  and 
•  is  a  term  that  accounts  for  both  the  (2J+l)-fold  Mj  degeneracy  and  the  Honl-London 

line'  strength  factor  for  the  specific  transition.35  The  ?i(l)  and  i?i(l)  data  are  fit  separately;  the 
results  of  these  fits  are  shown  in  Fig.  5.  •  The  fact  that  the  two  curves  do  not  coincide  is  due  to 
small  devajjons’  from  the  gas  phase,  Honl-London  line  strength  factors  in  the  matrix 
environment  From  the  two  fits  shown  in  Fig.  5  we  obtain  an  average  B  value  of  5  ±1  cm  , 
where  the  quoted  uncertainty  represents  the  difference  in  the  B  values  obtained. from  the  two  fits. 
This  value  is  consistent  with  the  rovibrational  assignment  presented  above.  Thus,  while  small 
deviations  from  gas  phase  behavior  are  observed,  the  appropriate  first-order  picture  that  emerges 
from  our  observations  is  one  of  nearly  free  rotating  DC1  in  the  pFk  matrix.  - 

-A-,- •.  One  more  prominent  deviation  from  the  gas  phase  spectrum,  however,  is  the  splitting  of 

.the  i?i(l)  transition  into  three  components.  While  the  origin  of  this  splitting  cannot  be 
unambiguously  determined  from  the  spectroscopic  data  alone,  we  speculate  that  it  arises  from  a 
lifting  .of  the  fivefold  44)  degeneracy  of  the  J=2  upper  rotational  state  that  is  accessed  in  this 
transition,  which  splits  this  state  into  three  Mj  components  (M/=0,  ±1,  ±2).4,5  This  splitting, 
which  may  signal  preferential  alignment  of  the  J=2  rotational  state  with  respect  to  the  pF^ 
crystal’s  c-axis,  presumably  arises  from  anisotropic  terms  in  the  orientational  potential  energy 
surface  experienced  by  the  J=1  DC1  state  in  the  pH2  matrix.  At  this  point^further  -theoretical 
modeling  is  needed  to  interpret  the  observed  .multiplet  structure,  but  the  first-order  picture 
remains  one  of  nearly  free  rotation  of  the  DC1  dopant.  Neither  the  Pi(l)  nor  the  (0)  transition- 
exhibits  any  visible  multiplet  structure.  This  is  interpreted  as  evidence  that,  to  within  the 
linewidth  of  these  transitions  (0.2  cm’1  FWHM),  the  pH)  matrix  does  not  lift  the  threefold 
degeneracy  of  the  J=\  state  in  either  vibrational  level. 

C.  HCl/pH2  rovibrational  assignment  .  )  . 

Based  on  the  results  for  DCl,  our  expectation  is  that  HC1  should  also  freely  rotate  in  solid 
pH2.  However,  due  to  HCl’s  larger  rotational  constant30  (5Hci=10.44  cm’1),  even  at  T=4-8  K  the 
ratio  of  the  J=\  to  J=  0  populations  should  be  only  1:175  (i.e.,  QI0t  ~  1).  Thus,  the  absorption 
spectrum  of  HC1  in  pH2  should  be  extremely  simple,  consisting  of  a  single  R\(0)  rovibrational 


transition  out  of  the  ground  J=  0  rotational  state.  Based  on  the  DC1  results,  this  transition  is 
expected  to  appear  approximately  10  cm'1  to  the  red  of  the  gas  phase  i?j(0)  transition.  , 
iiiP  thie-_surveyAScans  presented  in  Fig.  2,  a  strong  transition  in  vthe^HCl  doped  sample  ds 
evident  just  to  the  red  of  2900  cm'1.  As  noted  earlier,  this-  transition  displays  feflfl  isotopic- 
abundance  ratio  and  2.1  cm'1  splitting30  characteristic  of  natural  isotopic  abundance  HC1.  This 
feature^  is  assigned,  to  f?i(0)  based  on  its  absolute  transition  frequency  and  rpigt^'ely  broad 
lineshape.  The  Rx(0)  transition  for  H35C1  is  at  2894.2  cm'1,  redshifted  by  12.0  cm'1  from  the  gas 
phase  i?i(0)  transition30  at  2906.2464  cm'1.  If  this  transition  were  instead  assigned  to  the  Qx( 0) 
transition,  it  would  be  blueshifted  by  approximately  8  cm'1  from  the  gas  phase  transition  energy; 
this  seems  highly  unlikely  given  the  shape  of  the  pH2-HCl  intermolecular  PES  and-the  “size”  of 
a  single,  substitution  site.  In  addition,  Qi(0)  transitions  of  impurity  species  in  pH2.  are  usually 
very  sharp  because  they  involve  pure  vibrational  excitation  which  is  only  weakly  coupled  to  the 
pH2  solid.  In  contrast,  transitions  involving  J>0  states  are  generally  broader,  presumably  due  to 
them  -re  efficient  couplingbetween  the  low-energy  dopant  rotational  degrees  of  freedom  and  the 
phonon  modes  of  the  solid. 

,  v  ;  .  To  verify  this  rovibrational  assignment  and  to  better  characterize  the' purity  of  the  HC1 
rovibrational  quantum  states  in  pH2,  spectra  were  recorded  at  the  highest  possible  temperatures 
in  order  to  detect  transitions  from  thermally  populate^excited  HC1  rotational  states.  At  7=4.8  K 
the  ratio  of  the  gas  phase  integrated  intensities  of  i?,(l)  to  £,(0)  should  reach  1:263,  making  it 
just  possible  to  detect  fee  i?i(l)  transition  in  the  thickest,' most  concentrated  samples.  Indeed,  a 
broad  feature  is  observed  for  both  H35C1  and  H37C1  doped  pH2  solids  that  is  approximately  2 B 
higher  in  energy  than  fee  corresponding  Rx(0)  transition.  The  frequencies  and  linewidfes  of 

these  transitions  are  presented  in  Table  I. for  both  isotopomers. 

An  enlarged  region  of  fee  monomer  spectrum  for  fee  H35C1  doped  sample  in  fee  vicinity 
of  fee  transition  is  shown  in  Fig.  6.  The  different  traces  in  this  figure  show  fee  reversible 
temperature  dependence  of  fee  feature  assigned  as  ^(1).  The  i?i(l)  transition  is  most  evident  in 
trace  (ej,  which  is  fee  spectrum  recorded  at  T  «  10.  K  during  fee  destructive  sublimation  of  the 
sample.  Fitting  fee  temperature  dependence  of  fee  i?](l)  transition  intensity  to  Eq.  (1)  is  not 
possible  due  to  fee  negligibly  small  integrated  peak  intensities  in  most  spectra.  (The  temperature 
in  fee  subliming  sample  cannot  be  monitored  wife  precision  because  fee  solid  detaches  from  the 
substrate.)  However,  fee  temperature  dependence  of  this  peak  is  similar  to  feat  of  the  “hot”  DC1 
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transitions,  clearly  indicating  that  this  transition  originates  from  an  excited  state  at  a  significantly 
higher  energy.  Attempts  to  locate  the  corresponding  HC1  P\{\)  transition  were  unsuccessful 
because  the  more  concentrated  sample  (284  ppm  H35C1)  used  to  measure  these  weak  “hot” 
rovibrational  transitions  exhibits  significant  congestion  in  the  ?i(l)  region  from  (HC1)„  cluster 
peaks.  •• 

The  i?i(0)  and  Ri(l)  transitions  are  both  relatively  broad  under  these  conditions  but 
permit  the  average  rotational  constant  I?avg  of  HC1  in  solid  pH2  to  be  estimated.  For  both  HC1 
isotopoihers  the  calculated  5avg  is  9.1  cm  *,  which  is  about  7%  smaller  than  the  gas  phase 
values30  of  9.825  cm'1  for  H35C1  and  9.811  cm'1  for  H37C1.  Assuming  a  2 BAvg  rotational  spacing 
between  the  Ri(0)  transition  -and  the  vibrational  origin  for  this  band,  the  HC1  band  origin  v0  » 
Ri(0)  -  IBavg  =  2876.0  cm'1  is  redshifted  by  »  9.6  cm'1  from  the  gas  phase  value.  Thus,  the 
rovibrational  spectrum  of  HC1  in  pH2  closely  resembles  that  of  gas  phase  HC1,  except  that  the 
spectrum  is  redshifted  by  roughly  10  cm'1  and  exhibits  a  rotational,  spacing  that  is  .93%  of  the  gas 
phase  value. 

The  spectroscopic  data  presented  for  all  four  HC1  isotopes  are  therefore  consistent  with 
nearly  free  rotation  of  the  HC1  impurity  within  the  pH2  solid.. Thus,  the  gas  phase  rotational 
quantum  number  J  associated  with  the  total  angular  momentum  of  the  HC1  remains  an 
approximately  good  quantum  number  when  the  HC1  is  trapped  in  the  pH2  crystal.  The  nearly 
free  rotational  motion -of  the  HC1  is  presumably  a  direct  result  of  the  large  rotational  constant  of 
the  hydride  and  the  isotropic  solvation  site  that  the  pH2  presents  to  the  HC1  impurity. 

IV.  RESULTS  AND  ANALYSIS:  pH2  TRANSITIONS  INDUCED  BY  HC1  AND  DCI 
DOPANTS 

The  H2  <2i(0)  vibrational  transition  is  infrared  inactive  in  pure  solid  pH2  but  is 
observed13,23,36  in  pH2  solids  doped  with  low  concentrations  of  oH2.  The  mechanism  by  which 
oH2  dopants  induce  infrared  activity  associated  with  the  Qi(0)  transition  of  nearby  pH2 
molecules  has  been  studied  extensively;37  in  brief,  the  permanent  quadrupole  moment  of  the  oH2. 
impurity  generates  an  electrostatic  field  that  induces  transition  dipoles  in  pH2  molecules  in  the 
host  matrix.  This  quadrupole-induced-dipole  transition  is  inactive  in  pure  pH2  because  a  J=  0 
pH2  molecule  is  a  spherical  object  with  no  multipolar  electrostatic  field.  The  oH2 ‘induced  Qi(0) 
IR  absorption  transition  in  solid  pH2  doped  with  oH2  appears  at  4153.1  cm'1  and  exhibits  a 
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characteristic  lineshape  with  FWHM  0.4  cm'1  that  is  attributed  to  vibron  hopping  that  delocalizes 
the  v=l  excitation  throughout  the  hep  pH2  crystal.37  Because  this  transition  is  only  IR  active  if 
the  oH2  impurity  undergoes  a  change  in  its  Mj  quantum  number,  it  can  be  more  fully  described 
as  a Qj(0)  +  O0(  1)  transition;  this  notation  emphasizes  the  required  presence  of  the  oH2  dopant 
molecule  undergoing  the  Mj  orientational  transition.  The  oH2  dopants  themselves  also  undergo 
IITallowed  Oi(l)‘ transitions  visible  at  4146.6  cm'1. 

Strong  external  electric  fields  can  also  induce  £>i(0)  infrared  activity  in  solid  pH2,  in  this 
case  via  the  Condon  effect.17  However,  the  vibron  hopping  line-broadening  mechanism  is 
inactive  for  such  a  transition  because  the  inducing  field  is  external  to  the  pH2  crystal  lattice;' 
consequently  the  vibrational  excitation  remains  localized  on  an  individual  pH2  molecule  and  the 
transition  is  very  sharp:  it  appears  at  4149.70  cm'1  (with  FWHM  0.03  cm'1)  at  the  low-frequency 
end  of  the  vibron  band.  Similar  narrow  transitions  are  observed  in  solid  pH2  matrices  irradiated 
with  y-rays;18  these  transitions  are  thought  to  arise  from  a  Condon  effect  induction  mechanism  in 
which  the  inducing  freld  is  generated  by  a  large  number  of  immobile  charged  species  distributed 
throughout  the  pH2  solid.  To  the  extent  that  this  field  is  uniform  over  large  spatial  regions,  it 
mimics  the  effect  of  an  externally  applied  field. 

We  have  established  that  HC1  and  DC1  dopants  in  pH2  matrices  retain  a  good  rotational 
quantum  number  J,  and  that  the  vast  majority  of  these  dopants  are  in  the  .7=0  ground  rotational 
state  at  7—2.4  K.  Like  pH2  molecules,  HC1  and  DC1  molecules  with  J—  0  have  no  permanent 
-electrostatic  field;  consequently  these  dopants  cannot  induce  IR  activity  in  the  pH2  matrix 
through  mechanisms  based  on  electrostatic  induction.  However,  J=  0  impurities  can  generate  IR 
activity  in  the  pH2  matrix  via  an  alternative  mechanism  based  on  transition  dipole  moments 
induced  by  isotropic  exchange  and  dispersion  interactions,  which  are  active  even  for  7=0 
species.  ’  Consequently  it  may  be  possible  to  observe  6>i(0)  pH2  transitions  in  our  HCl  and 
DC1  doped  pH2  samples. 


A.  HC1  and  DC1  induced  Qi(0)  H2  transitions  in  solid  pH2  /vc(  (J) 

Survey  scans  of  the  H2  fundamental  region  (4100-4200  cm'1)  of  HCl^doped  pH2  solids 
are  shown  in  Fig.  7.  Traces  (a)  through  (d)  are  as-deposited  samples  recorded  at  temperatures 
near  2.4  K  for  pH2  samples  containing,  respectively,  only  residual  oH2,  HC1  dopants,  DC1 
dopants,  and  both  HC1  and  DC1  dopants.  The  very  weak  feature  just  discemable  in  trace  (a)  near 
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4153  cm'1  is  the  O\(0)  pH2  absorption  induced  by  residual  0H2  molecules  in  the  pH2  solid. 

The  broad  absorption  at  the  high  energy  end  of  the  spectrum  is  the  rising  edge  of  the  O^fOjH? 
phonon  sideband  which  involves  vibrational  excitation  as  well  as  phonon  excitation.  The 
extremely  sharp  features  near  4149.4  cm  1  in  traces  (b)  through  (d)  are  induced  by  the  HC1  and 
DC1  dopant  molecules.  The  frequencies  and  linewidths  of  these  features  are  summarized  in 
Table  II.  The  fact  that  these  features  are  observed  in  the  vicinity  of  the  H2  vibrational 
fundamental,  and  appear  at  nearly  the  same  frequency  for  both  the  HC1  and  DC1  doped -samples, 
supports  the  claim  that  they  involve  vibrational  excitation  localized  on  a  pHb  molecule  and  are 
not  associated  with  vibrational  or  rotational  excitation  of  the  impurity;  we  assign  these  features 
to  the  <9](0)  transitions  of  pH2  molecules  adjacent  to  the  HC1  or  DC1  impurity.  Traces  (c)  and  (d) 
also  show  the  i?2(0)  D35C1  arid  Dj7C1  overtone  rovibrational  transitions  that  serendipitously 
absorb  in  this  region. 

High-resolution  scans  of  the  induced  Q\(0)  pH2  feature  are  shown  in  Fig.  8, 
demonstrating  the  extremely  sharp  nature  of  this  feature  in  both  the  HC1  and  DC1  doped  solids. 
Traces  (a)  and  (b)  are  for  as-deposited  D35C1  and  Dj7C1  doped  pH2  samples,  respectively,  while 
traces  (c)  and  (d)  are  for  the  corresponding  as-deposited  HC1  isotopomer  doped  samples.  The 
linewidths  of  these  features  are  close  to  the  experimental  resolution  of  0.005  cm  allowing  for 
very  precise  determinations  of  the  transition  energies.  Note  that  given  the  narrowness  of  these 
features,  it  is  easy  to  resolve  the  0.07  cm"1  shift  between  the  HC1  and  DC1  induced  features.  The 
fact  that  these  features  depend  on  the  mass  of  the  dopant  species  is  consistent  with  our 
hypothesis  that  these  features  are  induced  by  the  dopants. 

The  features  shown  in  Fig.  8  appear  at  frequencies  very  close  to  the  Q\(0)  pH2  transition 
attributed  to  Coridon  effect  induction  in  solid  pH2  exposed  to  an  external  electric  field;  like  the 
Condon  effect  transitions,  these  features  are  also  very  narrow.  An  alternative  explanation  for  the 
features  depicted  irt  Fig.  8  might  therefore  involve  Condon-type  induction  by  the  superimposed 
dipolar  fields  of  the  (HC1)2  dimers  present  at  low  concentrations  in  the  as-deposited  samples. 
However,  annealing  of  these  samples  does  not  result  in  the  irreversible  growth  of  these  features. 
Instead,  additional  sharp  lines  at  nearby  frequencies  grow  in  upon  annealing;  these  lines  are 
attributed  to  induction  mechanisms  involving  HC1  dimers  and  larger  clusters.  Furthermore,  Fig. 
9  shows  that  in  as-deposited  samples,  the  integrated  intensity  of  the  narrow  features  shown  in 
Fig.  8  is  linearly  correlated  with  the  integrated  intensity  of  the  HC1  R\(0)  monomer  transition 


over  a  wide  range  of  HC1  concentrations.  In  contrast,  if  these  features  were  induced  by  (HC1)2 
dimers  (or  larger  clusters),  then  the  data  shown  in  Fig.  9  would  exhibit  nonlinear  character. 
Because  the  HC1  and  DC1  induced  features  shown  in  Figs.  7  and  8  are  present  in  as-deposited 
samples,  display  intensities  ,  linearly  correlated  with  the  intensities  of  known  monomer 
absorptions,  and  do  not  increase  in  intensity  upon  annealing,  we  conclude  that  they  are  induced 
by  HC1  or  DC1  monomers  and  not  by  clusters.  Because  nearly  all  of  the  HC1  and  DC1  monomers 
in  our  samples  are  in  the  J=  0  rotational  state,  we  assign  these  features  to  Qi(0)  H2  +  Q0(0)  HC1 
(or  DC1)  transitions;  here  the  notation  “+Q0( 0)”  emphasizes  the  fact  that  the  features  are  induced 
by  J-0  HC1  or  DC1  molecules. 

To  confirm  that  these  features  are  not  due  to  residual  oH2  molecules  in  the  solid,  samples 
of' HCl/pH2  were- intentionally  doped  with  higher -concentrations  of  oH2.  The  results  of 
intentionally  increasing  the  oH2  concentration  in  pH2  and  HCl/pH2  samples  are  shown  in  Fig.  10. 
Trace  (a)  shows  the  absorption  spectrum  at  7=2.4  K  of  an  as-deposited  33  ppm  H37Cl/pH2 
sample  also  containing  roughly  100  ppm  of  residual  oH2.  Trace  (b)  shows  the  spectrum  at  7=2.5 
K  of  an  as-deposited  pH2  sample  containing  2800  ppm  oH2  and  no  HC1;  this  spectrum  has  been 
rescaled  by  a  multiplicative  factor  of  0.5  for  purposes  of  comparison:  Trace  (c)  shows  the 
spectrum  at  7=2.4  K  of  an  as-deposited  pH2  sample  containing  both  52  ppm  H37C1  and  1200  ppm 
oH2. 

The  only  discemable  induced  feature  in  trace  (a)  is  the  6>i(0)  H2  +  Q0( 0)  HC1  induced 
absorption;  no  oH2  induced  H2  transitions  are  visible  in  this  spectrum.  In  trace  (b),  the  0\{\)  oH2 
pure  vibrational  transition  is  evident  at  4146.65  cm"1  and  the  oH2  induced  0i(O)  pH2  transition  is 
observed  at  4453.1  cm"1.  The  frequencies  and  linewidths  of  these  induced  H2  transitions  agree 
with  literature  values.13,23  Enrichment -with  oH2  [trade  (c)]  does  not  increase  the  intensity  of  the 
HC1  induced  feature  at  4149.4  cm-1,  thereby  ruling  out  the  possibility  that  this  feature  is  oH2 
dependent.  New  absorption  features  near  4139.65,  4143.11,  and  4149.66  cm'1,  marked  with 
asterisks,  appear  only  in  the  oH2  enriched  HC1  doped  sample;  assignment  of  these  features, 
which  are  presumably  due  to  oH2-HC1  clusters,  is  beyond  the  scope  of  the  present  paper. 

It  should  be  emphasized  that  the  J—  0  HC1  induced  (9i(0)  H2  transition  is  extremely  weak. 
The  main  reason  that  it  is  detected  in  the  current  experiment  is  that  the  transition  is  very  sharp, 
which  permits  it  to  be  detected  using  high-resolution  FTIR  methods.  In  fact,  if  the  DCl-induced 
transitions  were  only  a  factor  of  two  broader  (and  with  the  same  integrated  intensity)  they  would 
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probably  have  remained  undetected  at  the  current  signal-to-noise  levels.  The  very  narrow 
linewidth  of  the  Qi(0)  transition  is  intriguing;  as  we  discuss  in  more  detail  later,  it  suggests  that 
the  pH2  molecule(s)  involved  in  this  transition  are-decoupled  from  the  rest  of  the  pH2  solid. 

At  present  it  is  not  known  why,  as  Fig.  8  shows,  the  HC1  induced  Qi(0)  feature  is 
considerably  more  intense  than  the  DC1  induced  feature  (even  after  accounting  for  differences  in 
concentration  and  spectral  resolution).  The  origin  of  the  small  satellite  peak  to  the  blue  of  the 
Oi(0)  transition  in  the  HC1  doped  samples  is  also  not  clear  at  this  time. 

B.  Cooperative  pH2-HCI  transitions  in  solid  pH2 

In  both  the  HC1  and  DC1  doped  samples;  fairly  broad  features  are  observed  in  Fig.  7  at 
energies  slightly  above  that  of  the  Qi(0)  H2  transition.  The  frequencies  of  these  dopant  induced 
absorption  features  depend  on  the  rotational  constant  of  the  impurity.  For  example,  in  HC1 
doped  pH2,  the  first  broad  doublet  feature  is  approximately  20  cm*1  above  the  pH2  Qi(0) 
transition,  while  in  the  DC1  doped  sample  the  analogous  feature  is  only  10  cm*1  above  the  pH2 
Q\(0)  transition.  This  suggests  that  the  spacings  between  these  broad  features  and  the  induced 
Q\(0)  H2  peak  are  not  accidental,  but  rather  are  determined  by  the  J=0->1  rotational  spacing  of 
the  HC1  or  DC1  impurity.  Furthermore,  the  integrated  intensities  of  these  transitions  correlate 
linearly  with  the  intensity  of  die  Ri(0)  monomer  transitions  (see  Fig.  9).  Accordingly,  these 
transitions  are  assigned  to  Qi(0)  pH2  +  l?o(0)  HC1  (or  DC1)  cooperative  transitions  in  which  the 
dopant  molecule  undergoes  a  pure  rotational  transition  in  concert  with  the  pure  vibrational 
transition  of  the  pH2  molecule.  In  these  transitions,  a  single  infrared  photon  excites  both  the 
rotational  coordinate  .of  the  HC1  or  DC1  impurity  and  the  vibrational  coordinate  of  a  nearby  pH2 
molecule.  In  the  language  of  van  der  Waals  cluster  spectroscopy,  these  transitions  correspond  to 
combination  bands. 

The  detailed  ro vibrational  assignments  presented  in  Fig.  7  are  conducted  by  correlating 
the  gas  phase  rotational  spacing  of  the  dopant  with  the  frequency  difference  between  the 
cooperative  transition  and  the  impurity  induced  <2i(0)  peak.  For  example,  the  peak  labeled  <9i(0) 
H2  +  i?0(0)  HC1  in  trace  (b)  of  Fig.  7  is  approximately  25Hci  or  20  cm'1  higher  in  energy  than  the 
.  HC1  induced  <2i(0)  H2  transition.  Similarly,  the  first  peak  to  the  blue  of  the  DC1  induced  Q\(0) 
H2  transition  is  approximately  22?dci  or  10  cm*1  higher  in  energy  than  .  the  <2v(0)  transition. 
Absorption  features  blue  shifted  by  roughly  6B  from  the  Q\(0)  transition  are  also  observed  in 
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both  the  HC1  and  DC1  doped  pH2  samples;  Based  on  this  approximate  65  spacing  for  both 
isotopomers,  these  features  are  assigned  to  0\(Q)  H2  +  5o(0)  HC1  (DC1)  cooperative  transitions. 

In  the  DC1  doped  samples  the  Qi(0)  H2  +  50( 0)  DC1  transition  is  observed  near  4180  cm'1 
with  a  partially  resolved  multiplet  structure  that  is  presumably  due  to  anisotropic  interactions  that 
lift  the  Mj  degeneracy  of  the  DC1  J=  2  state.  In  the  case  of  HC1,  the  Oi(0)  H2  +  So(0)  HC1 
transition  appears  in  the  low-energy  portion  of  the  QR( 0)  phonon  sideband.  Thus,  in  order  to 
examine  the  lineshape  of  this  cooperative  transition,  we  recorded  the  absorption  spectra  of  two 
samples  with  the  same  amount  of  pH2  but  with  different  HC1  concentrations  and  subtracted  one 
spectrum  from  the  other.  Because  the  only  difference  between  the  two  samples  is  the 
concentration  of  HC1,  the  phonon  band  is  removed  by  subtraction  leaving  only  the  HC1 
induced-features,  as  shown  in  Fig.  1 1 .  This  figure  clearly  shows  the  Qi(0)  H2  +  S0(0)  H35Clpeak 
at  roughly  4209  cm'1.  Note  that  this  transition  is  somewhat  broader  than  the  corresponding  <3i(0) 
H2  +  S0( 0)  DC1  feature  shown  in  Fig.  7,  and  lacks  any  discemable  fine  structure.  This  is 
presumably  due  to  the  better  resonance  matching  of  the  HC1  rotation  with  the  pH2  phonon 
densityofstates^which  allows  for  faster  relaxation  of  the  HC1  rotation. 

.  The  Qi(0)  H2  +  i?o(0)  HC1  (DC1)  transitions  are  shown  on  an  expanded  scale  in  Fig.  12. 
The  difference  between  the  Qi( 0)  peak  and  the  average  of  the  frequencies  of  the  <2i  (0)  H2  +  R0( 0) 
HC1  doublet  feature  (weighted  by  the  approximate  2:1  intensity  ratio)  is  the  J=0-»1  rotational 


spacing  of  a  v=0  state  HC1  molecule  next  to  a  v=l  pH2  molecule.  This  gives  5=10.22  cm'1  and 
10.20  cm'1  for  the  H35C1  and  H37C1  isotopes,  respectively,  values  which  are  very  close  to  the 
Hj5Cl  gas  phase  value30  of  10.44  cm'1.  A  similar  analysis  gives  5=5.22  and  5.20  cm"1  for  D35C1 
and  D37C1,  respectively,  as  compared  with  the  corresponding  D35C1  gas  phase  value29  of  5=5.392 
cm'1.  These  5  values  compare  favorably  with  the  5Avg  values  determined  .from  the  direct 
infrared  spectroscopy  and  with  the  less  precise  value  of  5  determined  in  the  Boltzmann  analysis. 
However,  the  three  5  values  measure  slightly  different  quantities;  5Avg  is  averaged  over  HC1 

y-  5  .  " 

vibrational  states  the  5  value  determined  in  the  thermal  analysis  is  for  HC1  (v=0)  surrounded  by 
ground  vibrational  state  .  pH2,  and  the  5  value  determined  from  analysis  the  cooperative 
transitions  is  for  HC1  (v=0)  next  to  a  vibrationally  excited  pH2  molecule. .  Thus,  subtle 
differences  in  the  the  three  5  values  are  not  unexpected. 

The  spectra  shown  in  Fig.  12  are  for  the  isotopically  enriched  samples,  so" the  multiplet 
structure  observed  in  the  Q\(Q )  H2  +  R0( 0)  HC1  (DC1)  transition  is  not  due  to  different 
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isotopomers.  The  splitting  is  in  fact  a  direct  measurement  of  the  splitting  of  the  threefold, 
degeneracy  of  the  J=l  rotational  state  of  the  HC1  (DC1)  monomer  by  an  adjacent  vibrationally 
excited  pH2  molecule.  This  splitting  is  1.130  cm'1  and  1.156  cm’1  for  HC1  and  DC1,  respectively. 

'  No' splitting  of  the  J=1  rotational ’state  was  measured  for  HC1  surrounded  by  ground  ' vibrational 
state  pH2  molecules.  This  implies  that  excitation  of  a  neighboring  pH2  molecule  changes 
qualitatively  the  orientational  PES  of  the  HC1  (DCi)  impurity  in  the  matrix  environment.  The 
isotopic  dependence  of  the  observed  splitting  is  consistent  with  this  interpretation.  The  absolute 
magnitude  of  the  splitting  is  approximately  the  same  for  HC1  and  DCI,  which  means  that  relative 
to  the  J=1  to  J=0  free  rotational  spacing,  the  splitting  is  a  factor  of  two  larger  for  DCI.  A  larger 
relative  splitting  for  DCI  is  expected  because  the  smaller  rotational  constant  of  DCI  makes  it 
easier  to  perturb  the  dopant’s  rotational  motion  by  anisotropic  terms  in  the  PES.  We  will  return 
to  this  point  below. 

Another  striking  observation  is  the  difference  in  the  linewidths  for  the  Q\( 0)  pH2  pure 
vibrational' excitation  and  the  cooperative  transitions  Ci(0)  H2  +  Ro(0)  HC1  or  Q\(0)  H2  +  So(0) 
HC1.  The  linewidth  of  the  Qi(0)  +  R0(0)  cooperative  transitions  are  in  almost  quantitative 
•  agreement  with  the  linewidths  measured  for  the  corresponding  i?i(0)  transitions  of  the  HC1  and 
DCI  monomers.  This  again  is  interpreted  as  evidence  that  these  transitions  involve  excited 
rotational  states  whose  rotational  degrees  of  freedom  couple  more  strongly  with  the  phonon 
modes  of  the  solid  than  does  a  pure  vibrational  transition  The  fact  that  the  <3i(0)  H2  +  S'o(O)  HC1 
linewidths  are -broader  than  those  for  DCI  . is  because  the  cooperative  transitions  involving  HC1 
dopants  are  more  nearly  resonant  with  the  peak  of  the  pH2  phonon  modes,  resulting  in  faster 
dephasing  and  relaxation  of  HC1  (/=2)  compared  to  DCI  (J=2).  - 

Finally,  we  comment  briefly  on  the  relative  integrated  intensities  of  the  induced  Oi(0)  H2 
transitions  and  the  cooperative  Q\(0)  +  Ro(0)  and  Q\(0)  +  5o(0)  transitions.  As  Figs.  7  and  11 
show,  the  integrated  intensities  of  the  cooperative  transitions  are  much  larger  than  that  of  the 
pure  vibrational  Q\(0)  H2  transition.  This  is  because  two  different  induction  mechanisms  give 
rise  to  these  spectral  features.  As  we  have  explained,  the  IR  activity  of  the  pure  vibrational  O\(0) 
H2  transition  arises  from  transition  moments  induced  by  isotropic  overlap  and  dispersion 
interactions;  these  transition  moments  are  fairly  weak.38  In  the  cooperative  transitions,  however, 
the  upper  state  involves  a  J> 0  HC1  (DCI)  impurity  which  is  no  longer  isotropid  and^herefore^ 
■  exhibits  a  multipolar  electrostatic  field.  This  field  induces  a  vibrational  transition  dipole  in  the 
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neighboring  pH2  molecule  undergoing  the  0i(O)  transition,  so  that  the  cooperative  transitions  are 
much  stronger  than  the  pure  vibrational  transition. 

C.  HC1  perturbed  Si(0)_H2  transitions 

Another  induce^'  zero  phononpH2  transition  that  is  observed  in  these  studies  is  the  Si(0) 
H2  transition  near  4486.1  cm'1,  which  corresponds  to  rovibrational  excitation  of -a- -single' pHa 
molecule  (v=l<-0,  J=2<-0).  In  pure  pH2  solids  the  Si(0)  H2  transition  is  IR  active  because  the 
hep  crystal  lattice  lacks  inversion  symmetry.10'  Spectra  in  the  vicinity  of  the  Si (0)  transition  are 
shown  in  Fig.  13  for  both  HC1  and  DC1  doped  pHb  samples.  The  pure  pH2  transition  at  4486.1 
cm'1  is  nearly1  optically  black  for  these  3-mm-thick  pH2  samples,  but.  weak  transitions  are 
observed  slightly  to  the  red  of  this  feature  in  both  the  HC1  and  DC1  doped  samples.  The 
positions  of  these  features  do  not  change  to  within  experimental  uncertainty  for  the  individual 
35C1  and  37C1  isotopomers.  The  intensity  of  this  additional  feature  scales  linearly  with  HC1 
dopant  concentration.  Thus,  the  new  transition  is  assigned  to  the  £i(0)  transition  of  pH2  with  a 
neighboring  HCl  or  DC1  molecule. 

The  HCl  perturbed  *?i(0)  transition  is  shifted  5.4  cm'1  to  the  red  of  the  pure  pEk  Si(0) 
transition;  this  shift  arises  from  intermolecular  interactions  between  the  HC1  dopant  and  the  v=l, 
J- 2  state  of  H2.  .  Interestingly,  the  corresponding  redshift  for  the  DC1  perturbed  feature  is. 
significantly  larger  at  7.0  cm'1.  Such  a  large  isotope  effect  (1.6  cm'1)  for  a  transition  that  is 
localized  on  a  neighboring  pH2  molecule  is  unusual.  In  the  case  of  the  induced  <9i(0)  H2 
transition,  by  comparison,  the  spacing  between  the  features  induced  by  HC1  and  DC1  impurities 
was  only  0.07  cm'1.  The  large  isotopic  dependence  of  the  perturbed  Si(0)  pH2  transition  suggests 
that  the  J=  0  HC1  and  DC1  molecules  interact  differently  with  the  vibrationally  excited  J-. 2  pH2 
molecule.  Furthermore,  attempts  to  identify  cooperative  transitions  built  off  the  <Si(0)  transition 
that  involved  HC1  or  DC1  free  rotational  excitation  were  unsuccessful.  These  findings  are 
consistent  with  a  model  in  which  the  J-2  excited  rotational  state  of  the  pF^  molecule  strongly 
hinders  the  rotational  motion  of  the  neighboring  HCl  or  DC1  species.  If  this  is  the  case,  one 
would  expect  that  the  DC1  rotational  motion  would  be  more  strongly  perturbed  than  the  HCl 
rotational  motion  due  to  the  smaller  free  rotor  spacing  for  the  DC1  species.  Consequently  the 
shift  of  the  HCl  or  DC1  perturbed  Si(0)  transition  from  the  pure  pH2  iS'i(O)  transition  reflects 
changes  in  the  rotational  or  librational  zero  point  energy-  of  the  HCl  or  DC1  impurity.  Attempts 
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to  confirm  this  interpretation  by  modeling  the  rotational  motion  of  the  nominal  J—  2  H2  +  J- 0 
HC1  (DC1)  pair  are  currently  underway. 


IV.  DISCUSSION  : 

A.  HC1  and  DC1  dynamics  in  solid  pH: 

The  direct  infrared  excitation  spectra  of  HCf  and  DO  provide  quantitative  measurements 
of  how  the  vibration  and  rotational  dynamics  of  these  species  change  when  solvated  in  solid  pT^. 
For  HQ  and  DC1,  solvation  in  solid  pH2  redshifts  the, vibrational  frequencies  by  8.7  and  9.6  cm'1, 
respectively,  from  the  corresponding  gas  phase  values.  This  indicates  that  the  v=l  impurity  is 
-lowered  in  energy  upon  solvation  to  a  greater  extent  than  the.  v=0  impurity.  Such  an 
interpretation  agrees  with  the  IR  spectroscopy  of  pFt-HCl  van  der  Waals .  cluster,  which 
measured  a  -1.09  cm'1  shift  in  the  HC1  fundamental  vibration  upon  complexation.29  The 
redshifts  observed  here  therefore  provide  further  experimental  information  on  the  dependence  of 
the  pH2-HCl  potential  on  the  HQ  stretching  coordinate. 

The  spectra  also  provide  information  on  the  dependence  of  the  pH2-HCl  potential  on  the 
rotational  degrees  of  .freedom  of  the  HC1  monomer.  It  is  clear  from  the  spectroscopic  analysis 
presented  above  that  the  total  angular  momentum  quantum  number  of  the  HC1  and  DC1 
impurities  remains  reasonably  good  in  the  pff2  solid.  However,  the  J=0—>1  rotational  energy 
spacing  of  HQ  or  DC1  dopants  in  pH2  is  4-7  %  less  than  that  of  isolated  HC1  or  DC1  molecules  in 
the  gas  phase.  This  suggests  that  upon  solvation  in  pH^the  J=1  state  of  both  impurities  is 

y 

lowered  in  energy  by  a  greater  amount  than  the  J=0  state,  perhaps  because  the  J= T  state  can 
preferentially  sample  attractive  portions  of  the  PH2-HCI PES  shown  in  Fig.  1 . 

The  f?i(l)  transition  of  the  DC1  isotopomer  in  pH2  exhibits  both  a  reduced  rotational 
spacing  and  visible  multiplet  structure.  The  observed  splitting  is  interpreted  as  evidence  that  the 
J=2  rotational  state  accessed  in  this  transition  can  interact  with  the  pH2  matrix  through  higher 
order  anisotropic  terms  in  the  intermolecular  potential.  Thus,  while  the  J—  1  state  is  lowered  in 
energy  but  not  split  into  resolvable  sublevels,  the  J=2  state  is  both  lowered  and  .  split  into 
sublevels.  This  is  in  accord  with  elegant  symmetry  arguments  that  indicate  in  a  perfect  £>3h  hep 
crystal  site  the  magnitude  of  the  splitting  of  a  J-2  rotational  state  is  significantly  greater  thar^of  a 
J=  1  rotational  state.4,5  While  such  arguments  do  not  take  into  account  any  distortion  of  the  pH2 
crystal  structure  around  the  impurity,  they  provide  a  valid-  first-order  description  of  our 
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observations.  All  of  these  findings  point  to  a  recurring  theme  throughout  this  study,  which  is  that 
the  dynamics  of  the  solute  and  solvent  are  intimately  coupled. 

B.  Dynamics  of  pH2-HCl  pairs  in  solid  pH2 

The  cooperative  transitions  built  off  of  the  impurity  induced  C?i(0)  H2  pure  vibrational 
transition  provide  a  quantitative  energy  level  diagram  for  the  rotational  states  '  of  HCl  next  to  a 
single  vibrationally  excited  pH2  molecule  in  solid  pH2.  Figure  14  compares  the  gas  phase  HCl 
and  DC1  rotational  energy  levels  with  the  rotational  energies  determined  spectroscopically  from 
the  cooperative  transitions.  For  the  low  rotational  states  (7=1  and  J- 2)  of  HCl  and  DCL  the 
rotational  level  spacing  in  pH2  is  almost  identical  to  the  gas  phase.  This  again  validates  the 
assertion  that  the  HCl  and  DC1  species  retain  their  gas  phase  quantum  number  Jin  the  pH2  solid. 
Yet  there  is  also  a  clear  splitting  in  the  rotational  levels  in  the  condensed  phase  that  are  absent  in 
the  gas  phase,  which  is  thought  to  be  due  to  a  lifting  of  the  Mj  degeneracy  of  the  excited 
rotational  states  by  the  presence  of  a  neighboring  v=l,  J=0  pH2  molecule.  The  vibrational 
excitation  of  the  nearest  neighbor  pH2  molecule  breaks  the  symmetry  of  the  single  substitution 
site  and  defines  the  pair  axis  between  the  vibrationally  excited  pH2  and  HCl  as  the  new 
symmetry  axis.  We  can  develop  a  qualitative  picture  of  the  HCl  rotational  dynamics  under  the 
assumption  that  this  axis  is  the  Mj  quantization  axis;  we  will  use  the  label  k  to  represent  the 
projection  of  the  J=\  rotor  angular  momentum  onto  this  axis  in  order  to  emphasize  the  analogy 
between  the  rotational  dynamics  of  the  pH2-HCl  pair  in  solid  pH2  and  the  rotational  dynamics  of 
the  gas  phase  pH2-HCl  van  der  Waals  dimer.29 

In  the  case  of  the  J=1  state,  the  approximate  1 .3  cm"1  splitting  and  the  nearly  1 :2  intensity 
ratio  between  the  two  transitions  suggests  that  the  lower  energy  transition  (weaker)  accesses  the 
£=0  state  while  the  higher  energy  feature  (stronger)  accesses  the  £=±1  state.  The  same  energetic 
ordering  of  k  states  has  been  calculated  for  internal  rotor  states  of  the  gas  phase  pH2-HCl  dimer 
that  correlate  with  /=  1  HCl.40  In  the  gas  phase  dimer,  the  7hci=1,  k=  0  state  is  lower  in  energy 
because  the  -wavefunction  for  this  state  preferentially  samples  the  global  minimum  of  the  pH2- 
HC1  PES  in  which  the  HCl  points  toward  the  H2  ( i.e .,  the  global  minimum  shown  in  Fig.  1).  In 
contrast,  the  7hci=1,  k  =±1  internal  rotor  state  wavefunction  has  a  maximum  probability  density 
that  corresponds  to  the  HCl  bond  being  at  nearly  right  angles  with  respect  to  the  pair  axis, 
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corresponding  to  the  saddle  point  between  the  two  minima  in  Fig.  1;  thus  the  /hci-I,  £=±1 
internal,  rotor  state  has  a  higher  energy. 

If  this  model  is  correct,  the  J=  2  HC1  and  DC1  states  accessed  in  the  Q\(Q)  +  50( 0) 
cooperative  transitions  should  also  be  split  into  triplets  with  k=0,  ±1,  and  ±2.  However,  the 
splitting  of  the  J=  2  state  is  not  resolved  as  cleanly  as  that  for  the  J=l  state  due  to  broadening  of 
the  transition  and  dilution  of  the  oscillator  strength"  over  a  greater  number  of  sublevels.  The 
<9j(0)  +  S0(0)  transition  may  also  be  weaker  since  it  derives  the  bulk  of  its  intensity  from  the 
-  quadruple  moment  of  HC1  while  the  O\(0)  +  R0(0)  transition  depends  largely  on  the  dipole 
moment  of  HC1.  Nonetheless,  the  Oi(0)  +  S0(0)  shown  in  trace  (c)  of  Fig.  7  exhibits  two  mam 
peaks  that  are  separated  by  1.49  cm'1,  possibly  with  some  weaker  features  at  higher  energy.  In 
the  case  of  HC1  (Fig.  1 1)  the  corresponding  transition  is  rather  broad  (8  cm'1  FWHM),  which 
may  in  part  be  due  to  crystal  field  effects. 

In  addition,  since  the  impurities  are  confined  to  single  substitution  site^there  may  be  a 
strong  couplings  between  translational  and  rotational  degrees  of  freedom.  The  resultant  splitting 
of  the  HC1  Mj  rotational  sublevels  accessed  in  the  cooperative  transitions  will  therefore  be 
affected  by  rotation-translation  coupling  between  these  two  low  energy  motions.  Based  upon 
simple  particle-in-a-box  ideas,  the  “excited”  translational  states  should  "be  at  energies  comparable 
to  the  HC1  and  DC1  rotational  energy  spacings,  possibly  leading  to  strongly  coupled  motions. 

Due  to  the  extensive  zero-point  motion  of  the  pH2  solid,  a  rigorous  analysis  of  the 
coupling  between  the  three  translational  and  two  rotational  coordinates  of  the  HC1  dopant  must 
also  account  for  translational  motion  of  the  pH2  molecules  and  the  possible  correlations  between 
these  motions.  Such  a  rigorous  treatment  is  possible;  using  quantum  Monte  Carlo  (QMC) 
techniques^provided  an  accurate  intermolecular  PES  is  available:  The  data  presented  here 
describing  HC1  rotational  motion  in  the  pH2  matrix  provide  the  experimental  observables  that 
could  be  used  to  test  rigorous  QMC  calculations  using  a  trial  pH2-HCl  PES  that  includes 
dependence  on  the  intramolecular  stretching  coordinates.  Such  a  theoretical  treatment  of  the  HC1 
rotationkl  motion  in  a  single  substitution  site,  when  combined  with  the  high-resolution  data 
presented  here,  should  provide  a  full  microscopic  picture  of  the  quantum  dynamics  of  the  HC1 
impurities  in  solid  pH2. 

Another  important  finding  of  the  present  study  is  that  the  goodness  of  the  HC1 J  quantum 
number  is  dependent  on  the  J—  0  nature  of  the  pH2  solvent.  The  anomalous  isotopic  dependence 
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of  the  HC1  perturbed  H2  5i(0)  transition  (Fig.  13),  and  our  failure  to  observe  cooperative 
transitions  built  off  of  this  feature,  are  taken  as  evidence  that  the  J=2  H2  molecule  restricts  the 
rotational  motion  of  the  HC1  impurity,  transforming  it  from  one  of  free  rotation  to  libration. 

C.  /=  0  HC1  induced  Qi(0)  H2  transition 

" "  The  sharp  features  shown  in- Fig.  8  are  assigned  to  pure  vibrational  <2i(0)  excitation  of 
pH2  molecules.  Although  the  Ch(0)  transition  is  IR  inactive  in  pure  pFk,  it  becomes  active  in  the 
doped  pH:  solids  studied  here  through  the  induction -of  transition  dipole  moments  via  an  isotropic 
overlap  mechanism  which  is  the  only  induction  mechanism  operative  for  spherical  impurities. 
The  magnitudes  of  these  transition  moments  are  expected  to  be  particularly  sensitive  to  the  plTr- 
FTC1  distance;38  consequently  the  features  shown  in  Fig.  8  represent  a  rich  source  of  information 
about  the  microscopic  structure  of  the  HC1  and  DC1  substitution  sites  in  pHh  solids.  An  earlier 
report26  shows  that  spherical  atomic  impurities  such  as  Xe  also  induce  IR  absorption  near  4150 
cm'1,  presumably  arising  from  a  similar  mechanism.  It  is  interesting  to  note  that  analogous 
features  induced  by  J- 0  D2  or  HD  impurities  in  solid  pH2  have  never  been  reported,  despite 
numerous ’comprehensive  high-resolution  studies  of  isotopically  substituted  solid  hydrogens.13,23 
This  may  be  because  these  transitions  are  obscured  by  the  broad  £>i(0)  feature  induced  by 
residual  0H2  molecules  typically  present  in  these  samples. 

The  sharpness  and  absolute  frequency  of  the  gi(0)  transition  induced  by  HC1  or  DC1  may 
at  first  seem  surprising  because  the  Qi(0)  transition  induced  by  oH2  is  much  broader  (0.4  cm'1 
FWHM)  and  occurs  at  a  significantly  higher  frequency  (4153.10  cm-1).16,23  In  contrast,  the 
linewidth  of  the  HC1  induced  Qi(0)  feature  is  0.01  cm'1  and  substantially  redshifted  to  4149.39 
cm'1.  The  linewidth  and  redshift  of  this  transition  are  intimately  linked  because  as  the  vibrational 
frequency  of  the  pH2  molecule  next  to  the  impurity  becomes  detuned  from  the  vibrational 
frequency  of  the  other  pH2  molecules  in  the  matrix,  the  probability  for  vibron  hopping  to 
neighboring  molecules  drops  dramatically.  Thus,  once  the  Q\(0)  pH2  vibration  is  excited  it 
remains  strongly  localized  next  to  the  J=0  HC1  impurity  molecule,  which  constitutes  a  “vibron 
sink”  in  the  matrix. 


V.  CONCLUSIONS 

We  present  here  high-resolution.  IR  studies  of  solid  pEb  doped  with  HC1  and  DC1 
impurities.  .  By  utilizing  the  rapid  vapor  deposition  method  developed  by  Fajardo  and 
coworkers, 26,27  we  show  that  it  is  possible  to  synthesize  millimeter  thick,  optically  transparent 
samples  of  HC1  doped  pEb  solids  that  are  amenable  to  high-resolution  spectroscopic 
investigation."  This  approach  is  shown  to  produce  well-isolated  HC1  species  at  dopant 
concentrations  of  10  to  100  ppm.  Warming  the  sample  to  temperatures  above  4  K  permits 
diffusion  of  the  HC1  species' through  the  plL  solid  resulting  in  the  formation  of  (HC1}„  clusters. 
The  current  study  focuses  on  the  rotational  and  vibrational  dynamics  of  isolated  HC1  species  in 
the  pEb  solid. 

This  study  indicates,  as  expected,  that  1101  and  DC1  undergo  nearly  free  internal  rotation 
inside  the  pEb  solid.  This  is  taken  as  evidence  that  the  HC1  occupies  a  single  substitution  site  in 
the  hep  crystalline  structure,  because  such  a  site  is  a  nearly  isotropic  cavity  in  which  the  HC1 
molecule  can  freely  rotate.  For  both  HC1  and  DC1,  the  nearly  free  rotational  dynamics  are 
inferred  from  the  close  resemblance  between  the  condensed  phase  spectrum  and  the  hypothetical 
low  temperature  gas  phase  spectrum. 

In  addition;  a  weak,  extremely  narco  w^_we^k  absorption  feature  is  observed  near  4149.4 
cm'^which  is  assigned  to  a  pure  vibrational  Q\(0)  pEb  transition  induced  by  J=0  HC1  and  DC1 
impurites  through  an  isotropic  overlap  induction  mechanism.  Thus,  the  observed  frequency  and 
intensity  of  this  feature  is  sensitively  related  to  the  pEb-HCl  interaction  potential  and  to  the 
properties  of  the  HC1  doped  quantum  solid.  Interestingly,  similar  induced  transitions  have  been 
observed26  for  pH2  samples  doped  with  atoms  such  as  Xe,  but  never  before  been  reported  for  pEb 
solids  containing  J=0  molecular  impurities. 

The  present  assignment  of  the  pIL  <2i(0)  transition  is  confirmed  by  (i)  investigation  of 
0H2  enriched  samples  that  show  that  the  intensity  of  this  feature  is  independent  of  the  0H2 
concentration  and  (ii)  the  observation  of  cooperative  transitions  built  off  of  the  ptL  Q\(0).  These 
cooperative  transitions  involve  pure  vibrational  excitation  of  thepEb  molecule  in  concert  with 
pure  rotational  excitation  of  the  HC1  (or  DC1)  impurity,  and  exhibit  the  expected  dependence  on 
the  rotational  constant  of  the  impurity.  These  cooperative  transitions  are"  apparently  a  general 
phenomenon  associated  with  small  molecular  dopants  embedded  in  pEb  matrices,  as  we  have 
observed  them  for  pEb  solids  doped  with  a  wide  variety  of  impurities.41  Because  the  rapid  vapor 


deposition  method  used  here  makes  it  possible  to  incorporate  nearly  any  type  of  chemical 
impurity  into  solid  pH2,  these  cooperative  transitions  offer  new  insight  into  the  dynamics  of 
molecules  solvated  in  a  quantum  solid.  . _• 
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Table  I.  Transition  frequencies  for  HC1  and  DC!  rovibrationai  transitions  observed  in  solid  pH2 
at  T  «  2.4  K.  Linewidths  (FWHM)  are  given  in  parentheses...  All  values  are  in  units  of  cm'1 . 


transition 

H35C1 

"  H37C1 ' 

'  D35C1  "  ' 

D37C1 

^i(i) 

a 

a 

2072.183  (0.17) 

2069.273  (0.20) 

*i(0) 

2894.19(1.2) 

2892.06(1.1) 

2092.57  (0.16) 

2089.52  (0.17) 

2912.4(3.3) 

2910.3  (4.8) 

210131  (0.5) 
2102.59(0.6) 
2103.89(0.4) 

2098.29  (-— )b 
2099.60  (03) 
2100.98  (0.4) 

R-JO) 

a  _  i  . 

5663  .D  (2) 

i  i 

5659.9  (2) 

4121:05  (0.2) 

4115.17(0.3) 

a  obscured  by  cluster  absorptions 
b  overlapping  structure 


Table  II.  Transition' frequencies  for  HC1  and  DC1  induced  H2  rovibrationai  transitions  observed 
in  asrdeposited  solid  pH2  samples  at  T  »  2.4  K.  Linewidths  (FWHM)  are  given  in  parentheses. 
All  values  are  in  units  of  cm'1. 


Transition 

H35C1 

,  h37ci  •  -  .... 

D35C1 

d37ci  . 

Qi(0) 

4149.394(0.011) 

4149.391  (0.012) 

4149.318(0.02) 

4149.314(0.02) 

gi(0)  +  i?o(0){ 

4169.02  (0.21) 
4170.18  (0.70) 

4169.05  (--)a 
4170.16(0.9) 

4158.99(0.09) 

4160.16(0.18) 

4158.92,(0.15) 

4160.10(0.23) 

O1(0)  +  30(  0){ 

4209(8) 

‘  ;;  (™)b 

’  4179.77  (—  )a- 
4181.26  (— )a 

4179.71  (--f 
4181.28  ( — )a 

;  51(0) 

a  .  i  ^ 

4480.73  (0.27) 

4480.73  (0.27) 

4479.14(0.23) 

4479.13  (0.2) 

a  overlapping  structure 
b  not  determined 
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FIG.  1 .  Contour  plot  of  the  attractive  region  of  the  H2-HCI  ab  initio  intermolecular  potential 
with  contours  spaced  at  20  cm'1  intervals.  The  angular  orientation  of  the' Hi  is  assumed  to  be 
isotropic,  reducing  the  full  4D  intermolecular  potential  to  an  effective  2D  surface.  The  H2  and 
HC1  bond  lengths  are  held  fixed.  Note  the  nearly  circular  zero  contour.  * 


wavenumber  (cm’1) 


FIG.  2.  HCl/pH2  and  DCl/pHb  survey  IR  absorption  spectra  with  Cl  present  at  natural  isotopic 
abundance.  The  upper  panel  shows  spectra  (0.0075  cm'1  resolution)  of  a  2.9-mm-thick  88  ppm 
HCl/pH2  solid  at  T  —  2.4  K;  trace  (a)  is  for  the  as-deposited  sample,  while  brace  (b)  is  for  the 
same  sample  recooled  to  T=  2.4  K  after  annealing  at  T  =  4.8  K  for  1  hour.  The  lower  panel 
shows  spectra  (0.01  cm'1  resolution)  of  a  3.0-mm-thick  102  ppm  DCl/pH2  solid  at  T-2.5  K; 
trace  (a)  is  as-deposited,  while  trace  (b)  is  taken  at  T  =  2.4  K  after  annealing  at  T  =  4.9  K  for  1 
hour.  In  both  samples  the  annealing  process  results  in  irreversible  growth  of  the  cluster 
absorptions. 
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FIG.  3.  Temperature  dependence  of  D35Cl/pH2  monomer  fundamental  absorptions  (0.1  cm'1 
resolution).  The  sample  is  a  3.1-mm-thick  49  ppm  D35Cl/pH2  solid;  trace_(a)  is  for  the  as- 
deposited  sample  at  T  =  2.5  K,  trace  (b)  is  during  annealing  at  T=  4.9  K,  trace  (c)  is  for  the 
sample  retooled  to  T=  2.5  K,  and  (d)  is  for  the  sample  warmed  up  to  T  =  3 .6  K.  The  (DCl)n 
cluster  absorptions  . grow  irreversibly  in  intensity  upon  initial  annealing,  but  the  ?i(l)  and  i?i(l) 
DC1  monomer  features  show  a  reversible  temperature  dependence.  v 


FIG.  4.  Temperature  dependence  of  D37Cl/pH2  monomer  fundamental  absorptions  (0.1  cm'1 
resolution).  The  sample  is  a  3.0-mm-thick  56  ppm  D37Cl/pH2  solid;  trace  (a)  is  for  the  as- 
deposited  sample  at  T=  2.6  K,  trace  (b)  is  during  annealing  at  T=  4.8  K,  trace  (c)  is  for  the" 
sample  recooled  to  T =  2.6  K,  and  trace  (d)  is  for  the  sample  warmed  up  to  T=  3.6  K. 


FIG.  5.  Temperature  dependence  of  the  integrated  intensity  of  the  Pi(l)  (circles)  and  i?i(l) 
(squares)  DC1  rovibrational  transitions  in  DCl/pEb  samples.  The  solid  lines  represent  separate 
fits  of  the  Pi(l)  and  J?i(l)  data  to  Eq.  1;  see  the  text  for 'further  details. 
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FIG.  6.  H35Cl/pH2  monomer  fundamental  absorptions.  Spectra  for  two  different  samples  are  ' 
shown  (1.0  cm"1  resolution).  Trace  (a)  is  for  an  as-deposited  3.0-mm-thick  30  ppm  H35Cl/pH2 
solid  at  T=  2.4  K;  the  spectrum  has  been  rescaled  by  a  multiplicative  factor  of  0.2.  Traces  (b) 
through  (e)  are  for  a  3.0-mm-thick  284  ppm  H35Cl/pH2  solid;  trace  (b)  is  for  the  as-deposited 
sample  at  T=  2.4  K,  trace  (c)  is  during  annealing  at  T—  4.8  K,  trace  (d)  is  for  the  sample 
recooled  to  T  =  2.4  K,  and  trace  (e)  is  during  the  destructive  sublimation  of  the  sample  at  T  ^  10 
K. 

•s. 
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FIG.  7.  Survey  of  HC1-  and  DCl-induced  IR  activity  in  solid  pH2.  Trace  (a)  is  the  spectrum  (0.1 
cm'1  resolution)  of  a  3 . 1  -mm-thick  as-deposited  120  ppm  oFI2/pFI2  solid  at  T  =  2.4  K.  The  very 
weak  feature  near  4153  cm"1  is  the  Oi(0)  H2  absorption  induced  by  residual  oH2  molecules,  the 
broad  feature  at  the  right  side  of  the  spectrum  is  the  rising  edge  of  the  Or( 0)  H2  phonon  sideband. 
Trace  (b)  is  for  the  as-deposited  88  ppm  HCl/pH2  sample  depicted  in  the  upper  panel  of  Fig.  2. 
Trace  (c)  is  for  the  as-deposited  102  ppm  DCl/pH2  sample  depicted  in  the  lower  panel  of  Fig.  2. 
Trace  (d)  is  for  an  as-deposited  2.9-mm-thick  pH2  sample  doped  with  50  ppm  HC1  and  100  ppm 
DC1  at  T  =  2.5  K;  0.01  cm'1  resolution.  The  sharp  features  in  traces  (b)  through  (d)  near  4149.4 

cm'1  are  the  <3i(0)  H2  absorption  transitions  induced  by  the  presence  of  HC1  and  DC1  dopant 

$ 

molecules.  The  labels  “+  R0( 0)  DC1,”  “+  R0( 0)  HCF’and  “+  50(0)  DC1”  are  shorthand  for  the 
2i(0)  H2  +R0(0)  DC1,  2i(0)  H2  +  i?o(0)  HCl,  and  Q{(0)  H2  +  5o(0)  DCl  cooperative  absorptions 
in  which  the  dopant  molecules  undergo  a  pure  rotational  transition.  Traces  (c)  and  (d)  also  show 
the  i?2(0)  DC1  vibrational  overtone  transition. 


FIG.  8.  Isotope  dependence  of  the  Q}(0)  H2 1  Qo(0)  HOI  induced  IR  absorptions.  Trace  (a)  is 
the  spectrum  (0.01  cm'1  resolution)  of  an  as-deposited  3.0-mm-thick  49  ppm  D35Cl/pH2  solid  at  T 
=  2.5  K.  Trace  (b)  is  the  spectrum  (0.01  cm'1  resolution)  of  an  as-deposited  3.0-mm-thick  56 
ppm  D37Cl/pH2  solid  at  T=  2.5  K.  Trace  (c)  is  the  spectrum  (0.005  cm’1  resolution)  of  an  as- 
deposited  3.0-mm-thick  90  ppm  H35Cl/pH2  solid  at  T  =  2.4  K.  Trace  (d)  is  the  spectrum  (0.005 
cm'1  resolution)  of  an  as-deposited  2.9-mm-thick  94  ppmH37Cl/pH2  solid  at  T=  2.4  K. 
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FIG.  9.  Correlation  plot  of  integrated  intensities  of  direct  and  induced  ER.  absorptions  in  dilute, 
as-deposited  HCl/pF^  solids.  The  circles  are  for  the  Q\(0)  H2  +  Ro(0)  HCI  transition  near  4170 
cm'1,  the  triangles  are  for  the  sharp  <9i(0)  H2  transition  near  4149.4  cm'1;  the  latter  have  been 
rescaled  by  a  multiplicative  factor  of  10.  • 
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FIG.  10.  0H2  and  7C1  induced  IR  activity  in  solid  pH2.  Trace  (a)  is  the  spectrum  (0.005  cm'1 
resolution)  of  an  as-deposited  3.0-mm-thick  33  ppm  H37Cl/pH2  solid  at  T—  2.4  K.  Trace  (b)  is 
the  spectrum  (0.008  cm'1  resolution)  of  an  as-deposited  4.0-mm-thick  2800  ppm  0H2/PH2  solid  at 
T  =  2.5  R;  this  spectrum  has  been  rescaled  by  a  multiplicative  factor  of  0.5.  Trace  (c)  is  the 
spectrum  at  T  =  2.4  K  (0.01  cm'1  resolution)  of  an  as-deposited  3.0-mm-thick  pH2  solid  doped  ’ 
with  52  ppm  ppm  H37C1  and  1200  ppm  0H2.  The  features  marked  by  asterisks  near  4139.65, 
4143.1 1,  and  4149.66  cm'1  only  appear  in  the  sample  doped  with  both  HC1  and  0H2. 


FIG.  11.  H35Cl-induced  ER.  activity  in  solid  pH2  (0.005  cm'1  resolution).  Trace  (a)  depicts  an  as- 
deposited  3. 0-mm-thick  30  ppm  H35Cl/pH2  solid  at  7  =  2.4  K;  trace  (b)  depicts  an  as-deposited 
3.0-mm-thick  284  ppm  H35Cl/pH2  solid  as-deposited  at  7=  2.4  K.  The  lower  panel  shows  the 
direct  subtraction  of  trace  (b)  from  trace  (a).  Cancellation  of  the  broad  Qr(0)  phonon  sideband 
reveals  the  weak  Q i(0)  H2  +  5o(0)  HC1  cooperative  absorption  near  4209  cm'1. 
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FIG.  12;  Isotope  dependence  of  the  Oi(O)  H2  +  Ro(0)  HC1  (DC1)  induced  IR  absorptions.  Trace 
(a)  is  the  spectrum  (0.01  cm'1  -resolution)  of  an  as-deposited  3.0-mm-thick  49  ppm  D35Cl/pH2 
solid  at  T- 2.5  K.  Trace  (b)  is  the  spectrum  (0.01  cm'1  resolution)  of  an  as-deposited  3.0-mm- 
thick  56  ppm  D37Cl/pH2  solid  at  T=  2.5  K.  Trace  (c)  is  the  spectrum  (0.005- cm'1  resolution)  of 
an  as-deposited  3.0-mm-thick  90  ppm  H35Cl/pH2  solid  at  T  =2.4  K.  Trace  (d)  is  the  spectrum 
(0.005  cm'1  resolution)  of  an  as-deposited  2.9-mm-thick  94  ppm  H37Cl/pH2  solid  a \T=  2.4  K. 


FIG.  13.  Isotope  dependence  of  the  Si(0)  H 2  +  Qo(  0)  HC1  (DC1)  induced  IR  absorptions.  Trace 
(a)  is  for  the  49  ppm  D35Cl/pH2  sample  depicted  above  in  Fig.  12  (a).  Trace  (b)  is  for  the  56  ppm 
D37Cl/pH2  sample  depicted  in  Fig.  12  (b).  Trace  (c)  is  for  the  102  ppm  DCl/pH2  sample  depicted 
in  trace  (a)  of  the  lower  panel  of  Fig.  2-.- Trace  (d)  is  for  the  90  ppm  H35Cl/pH2  sample  depicted 
in  Fig.  12  (c).  Trace  (e)  is  for  the  94  ppm  H37Cl/pH2  sample  depicted  in  Fig.  12  (d).  Trace  (f)  is 
for  the  88  ppm  HCl/pH2  sample  depicted  in  trace  (a)  of  the  upper  panel  of  Fig.  2. 


70  i 


60  H 


E  50 
>* 

B  40 
CD 
c 
CD 

ffl'  30 
C 
.2 

J2  20 
o 


10  A 


0  J 


J*. 2  ___ 

FWHM 


J=1 


J=.Q' 


T 


gas  pH2 

HCI 


:J=. 2 


J=  1 


J=0 


gas  pH2 

DC  I 


Fig.  14.  Rotational  energy  levels  of  HCI  and  DC1  in  the  gas  phase  (gas)  and  solvated  in  solid 
pH2  with' one  vibrationally  excited  nearest  neighbor  (jaFb):'  Both  solvated  species  display  an 
approximately  1  cm*1  splitting  of  the  J=  1  rotational  energy  level.  The  J- 2  rotational  state  is  split 
into  at  least  two  resolvable  sublevels  for  DC1,  but  only  a  broad  energy  range  can  be.  specified  for 
HCI,  specified  by  the  FWHM  of  the  observed  spectral  feature. 


